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The effect of focusing and defocusing nonlinearities on Anderson localization in highly nonlocal 
media is theoretically and numerically investigated. A perturbative approach is developed to solve 
the nonlocal nonlinear Schroedinger equation in the presence of a random potential, showing that 
nonlocality stabilizes Anderson states. 



Disorder and nonlinearity are two leading mechanisms 
promoting wave localization. On one hand, a sufficient 
strength of disorder fosters the transition from a diffusive 
regime to a wave-function exponentially decaying over a 
characteristic distance Z; a scenario commonly referred 
to as Anderson localization [l|-[lo|. On the other hand, 
in a nonlinear medium, diffraction, or dispersion, can be 
compensated by the nonlinear properties; when the char- 
acteristic lengths of these phenomena (nonlinearity and 
diffraction) are comparable, localized solitary waves, or 
solitons, emerge. It is well known that nonlocality in the 
nonlinear response can largely affect the localized wave- 
forms, depending on the degree of nonlocality a [iMi. 
However the interplay between the disorder induced lo- 
calization length I and the characteristic length of non- 
locality a has never been considered before. Here we 
analyze the effect of a nonlocal nonlinearity on Ander- 
son localization; we show that in the framework of the 
highly nonlocal approximation, it is possible to derive 
closed form expressions to describe the role of nonlin- 
earity on Anderson localizations, and that these states 
become more stable when the degree of nonlocality in- 
creases, meaning that the power needed to destabilize 
them increases with a/ 1; this result unveils a fundamen- 
tal connection between nonlocality and disorder. 
Model — The nonlocal nonlinear Schroedinger equation 
reads as 

/ + CO 
X{x' - x)\'ip{x')\^dx' (1) 
-oo 

where ^ — 'ip[x,t)^ ^{x) is a random potential and 
x{x) is the response function of the nonlocal medium 
normalized such that J x{x)dx = 1; s ± 1 corresponds 
to a focusing (s = 1) or defocusing (s = —1) nonlin- 
earity. Eq. ([T]) applies to a variety of physical prob- 
lems, including nonlinear optics and Bose Einstein con- 
densation We define the "unperturbed" Hamilto- 
nian as — —d^ + V{x); its eigenstates are written 

as HoTpn = (3ni^n with {^n,^m) = ^nra- Hq SUStaiuS 

exponentially localized states, corresponding to negative 
eigenvalues /3„. The fundamental state can be approx- 
imated by ■0o(a^) = -^^e~^^~^°^/'^ , where the average lo- 
calization length I is determined by the strength of the 
random potential Vq , and xq is the location of the eigen- 
function with eigenvalue Pq. In the following, without 



loss of generality, we assume that the horizontal axis has 
been shifted such that xq = 0. From Hq, one roughly 
finds the link between the eigenvalue and the average lo- 
calization lenght / by a taking the average wavefunction 
V'o: Wo) = {{ipO:Hoipo)) — —1/^^; the lowest energy state 
has the highest degree of localization. 
Highly Nonlocal Limit — The nonlocality is described by 
x(x), which is typically bell-shaped with a characteris- 
tic length a. For an average localization length I much 
shorter than the nonlocality degree cr of the medium, the 
response function x{^' ~ can be expanded around the 
localization center, x' — 0, and one has in ([1]) 

/ + 00 
X[x' -x)\i,[x')\^dx' = Px{x). (2) 
-oo 

In this highly nonlocal limit (HNL), the nonlinearity 
can be treated as an interaction Hamiltonian Hint = 
— sx(a;)P, where P — /J*^ \i\){xy\'dx is overal energy, 
or beam power, which is conserved during evolution af- 
ter Eq.([T]). This limit is valid in the regime of a wave 
dominated by a single localization, which, without loss 
of generality, is taken centered at = 0. This also holds 
true as far as during the dynamics, additional localiza- 
tion are generated among those located in proximity of 
X = 0. Conversely, if two distant localizations are ex- 
cited, x(^) will be composed by two nonlocal responses 
centered in the two Anderson states, this case will be 
investigated elsewhere. The HNL allows to apply the 
standard perturbation theory of quantum mechanics for 
deriving closed form expressions for the effect of non- 
linearity on the Anderson states. We write the field as 
an expansion in P, V = ^J~P (V'o + -PV'^"^^ + -P^V-"^^^ + ••■)' 
where we take at the leading order i)) = \/P^o to focus 
on the effect of nonlocal nonlinearity on the fundamental 
state. We obtain the correction to the Anderson ground 
state (at second order in P) eigenvalue: 

/3o(P) = /3o-.Pxoo + P^E#^' (3) 

with the matrix elements of the nonlocality given by 
Xnm — J xix)ipn{x)Tpm(x)dx . As the degree of nonlocal- 
ity increases, x(^) can be treated as a constant in Xnm, 
such that Xmn — x{^)^nm- This shows that the pertur- 
bation to the Hamiltonian is diagonal, therefore, in the 
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HNL, the effect of the nonlocahty is to shift the eigen- 
value such that (3o{P) = 13q — sx{0)P, where x(0) depends 
on the specific x{x)j exphcitly contains the degree of 
nonlocahty. The overall Hamiltonian H — Hq + Hi^t is 
diagonal in the same states of Hq, hence Anderson local- 
izations turn out to be eigenstates also in the presence of 
the nonlocahty. These localizations are hence expected 
to be extremely robust with respect to the nonlinearity. 
Instability of Anderson states — The effect of nonlinear- 
ity on Anderson states becomes relevant when the term 
linear in P is comparable with /3o (higher order correc- 
tions vanish in the HNL); this allows to define through 
Eq. ([3]), the critical power Pc : 
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(4) 



For a defocusing medium (s ~ ^1), this is the power 
needed to change the sign of the eigenvalue l3o{P), from 
negative to positive; such that the localization is de- 
stroyed. Conversely, for a focusing medium (s = 1), this 
can be interpreted as the power were the average degree 
of localization is strongly affected by nonlinearity, indeed 
as 1{P) ^ l/y/{l3{P)), one has 



1{P) = 



v/l + sP/(Pc) 



(5) 



such that at critical power, the localization length is re- 
duced by a factor \/2 for the focusing case s — 1, and 
diverges for the defocusing case s = — 1. This shows that 
for P > Pc no localized states are expected for s = — 1 , as 
the corresponding eigenvalue changes sign. This trends 
applies as far as additional effects, like the excitation of 
further localizations, occur. From Eq. ([4]), we obtain the 
expression for P^ in the HNL, P^ = A = |/3o|/|x(0)| 
with (Pc) = a/2P [x(0) = l/2a for an exponential non- 
locality] . Because of the Cauchy-Schwarz inequality, one 
has \4>o{x)\'^xix)dx < 1, and one readily sees that 

for a finite nonlocahty Pc < Pc', as the nonlocahty in- 
creases the power needed to destabilize the Anderson 
states grows. 

A useful measure to quantify the effect of a nonlocal non- 
linearity on Anderson states is the residual value of /?(P) 
at the critical power Pc, which can be written as 



XooP(l-/3„//3o) 



(6) 



As it is determined by the off-diagonal elements Xno, 5 13 
can be taken as the "residual coupling" due to the non- 
local nonlinearity, which vanishes in the HNL limit. In 
Figlll we numerically show [for an exponential x{^)\ that 
5(3 goes to zero when increasing a: as the nonlocahty in- 
creases the nonlinear coupling of V'o with other states is 
moderated, hence it tends to behave as an eigenstate of 
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FIG. 1: (Color online) Left axis: critical power for ex- 
ponential [Pc = |/?o|(2(T -I- I), dashed-dotted line), Gaus- 
sian {Pc = ii/3o| exp(— (T^//^)/i5r/c(cr//), continuous line), 
quadratic (Pc = |/^o|/|x(0)|(l — ^^/2(t'^), dots) and rectangu- 
lar (Pc = 2cr|^o|/(l - e"^'"/'), diamonds) response functions, 
for a single realization of the disorder; right, residual cou- 
pling (5/3 Vs degree of nonlocality, averaged over 10 disorder 
realizations (Vo = 1). 



the system even if nonlinearity is present (however, its 
degree of localization may be largely affected) . 
Nonlocal Responses — We analyzed a few specific case 
of response function x(a;) (rectangular x = l/(2o') 
for |a;| < cr, x = elsewhere; exponential x{^) — 
g-|:r|/<T/(2cr); Gaussian x{x) = e'^'/'^'/vw^; and 
quadratic xi^) = x(0) + Xaa;^). In Fig. [U we report the 
behavior of the critical power as a function of the non- 
locality degree a for the analyzed response functions. In 
all of these cases, the calculated critical power is linearly 
dependent on the unperturbed eigenvalue of the state. 
So, the higher the strength of the disorder Vq, the lower 
/3o, the higher the power needed to affect the localiza- 
tion. Furthermore, we emphasize that the critical power, 
depending on the Anderson eigenvalue, has a statistical 
distribution depending on the disorder configuration. 
Numerical Results — We numerically solved Eq.([T]) for 
fixed disorder configurations: we consider a Gaussianly 
distributed V{x) with zero mean and standard deviation 
Vo, first obtain the eigenstates, then by using a pseudo- 
spectral Runge-Kutta algorithm, we evolve the Ander- 
son localizations in a nonlocal medium with a given x{^) 
[an exponential response hereafter, similar results are ob- 
tained for other xi^)]- Fig. [2] shows the dynamics of the 
ground-state intensity: for P < Pc the state remains al- 
most unperturbed (Fig. [2l^a), attractive case; Fig. [2jc), 
repulsive case). By increasing the power beyond the 
critical threshold, we observe two different phenomena. 
In the focusing case, the state becomes more localized 
Fig-IHb). At higher powers, a temporal beating pattern 
can be observed. This is mainly due to the coupling of 
the Anderson ground state with other localizations. In 
fact, for cr = 10 the highly nonlocal limit is not satisfied, 
and a residual nonlinear coupling between the Anderson 
ground state and additional localized modes is present 
(see Fig. [H 2nd y-axis), and causes the observed oscil- 
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new experimental investigations, can be extended to sev- 
eral related problems, as quantum phase diffusion and 
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FIG. 2: (Color online) Evolution of the ground-state intensity 
for a fixed disorder realization, for a = 10, Vo = 10; focusing 
(attractive) case s = 1, for P = 0.04Pc (a) and for P = APc 
(b); defocusing (repulsive) case s = — 1, for P = 0.04Pc (c) 
and for P = 4Pc (d). The superimposed dashed hue in the 
panels represents the ground-state at the initial time. 
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FIG. 3: (Color online) Wavefunction waist for focusing (a), 
and defocusing (b) cases, for P = 0.04Pc (dash-dotted line), 
P = OMPc (continuous line), P = 2Pc (dotted line), P = 4Pc 
(dashed line), (cr = 10, Vq = 10). Results averaged over 10 
disorder reahzations. 



lations. In the defocusing case, we observe the break- 
ing of the Anderson localization, as expected (Fig. [2l^d)). 
Fig. [3] shows the localization length of the ground state 
for various powers. In the focusing case, the localization 
length decreases and the beating pattern is observed. For 
s = — 1, the perturbation delocalizes the eigenmode as 

P>Pc. 

Conclusions — We reported on a theoretical analysis of 
the effect of a nonlocal nonlinearity on disorder induced 
localization. We derive explicit formulas to predict the 
critical power for destabilizing the Anderson states, in 
quantitative agreement with numerical simulations. We 
have shown that an increasing degree of nonlocality pro- 
duces a substantial growth of the power needed to desta- 
bilize Anderson states, which turn out to be very robust 
with respect to nonlinear effects. This can also be ex- 
plained by the fact that nonlocality reduces the coupling 
between Anderson states. These results may stimulate 



coherence [17|, ultrashort pulses in fibers [18[, second 
harmonic generation [l^ and Bose-Einstein condensation 



and furnish a novel theoretical framework for the 
interplay of nonlinearity and disorder (20j . 
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